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Abstract. - We analyse Bloch oscillations (i.e., oscillations induced by a static force) of 
strongly interacting Bose atoms beyond the hard-core bosons model. It is shown that residual 
interactions between the atoms modulate Bloch oscillations, where the type of modulations 
depends on the magnitude of the static force. 



1. Recently much attention has been paid to the transport of cold atoms in optical lattices 
[1-12]. One of these problems is Bloch oscillations (BO) of atoms in a static, for example, 
gravitational field (see review [6] and references therein). Besides their own interest, these 
studies pursue the aim of realizing (and understanding) the ordinary and super-conductivity 
of cold atoms [8, 12]. 

The present work is devoted to BO of Bose atoms in deep anisotropic optical lattices, 
where the system can be described by the ID Bose-Hubbard model 

H = H + Fd^lhi , (1) 
l 

H = -^Y1 + h - c ) +yE "'("J - !) • ( 2 ) 

i ( 

In Eqs. IlllL'll a[ and di are the bosonic creation and annihilation operators, hi = d\di the 
number operator, J the hopping matrix element, W the on-site interaction energy, d the 
lattice period, and F the magnitude of static force. For vanishing interactions the static 
force in (JIJ induces coherent, perfectly periodic oscillations of the atoms with the Bloch 
frequency lub — Fd/Ti. Then the main question is of how non-zero interactions do affect these 
oscillations. 

For a weak to moderate interaction constant, W < J, BO of Bose atoms in ID lattices were 
analysed in Ref. [4]. It was found that atom-atom interactions destroy the coherence in the 
system and BO typically decay within a few Bloch periods. (An exception is the case of strong 
forcing, Fd 3> J, where BO of interacting atoms are quasiperiodic [5].) One might naively 
expect that stronger interactions would cause an even faster decay of BO. This is, however, 
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not the case because for W 3> J the system approaches the hard-core bosons (HC-bosons) 
limit, 

H H c = - 2 + h - c ) ' ( 3 ) 

where the additional constraints (a;) 2 = {a\) 2 = prohibit double occupancy of a single well. 
The system (0 is completely integrable and its Hamiltonian can be diagonalised by using the 
canonical transformation bk = (1/VL) J2i exp(i27rfc^/L)a;, i.e., by going from the Wannier to 
Bloch basis. Then it is easy to show that BO of HC-bosons do not decay [see Eq. (0 below]. 

The above result refers to the formal limit W — > oo. In practice, however, one deals with 
a finite interaction constant, which implies the presence of residual interactions between Bose 
atoms. It is the primary aim of this work to study the effect of these residual interactions on 
the Bloch dynamics of strongly interacting Bose atoms. 

2. We begin with the effective Hamiltonian for strongly interacting Bose atoms( 1 ) . As- 
suming periodic boundary conditions the Hilbert space Hhc of the system © is spanned by 
the Fock states |n) = \rii, . . . ,til) with the occupation numbers n; being ether zero or unity. 
To take the residual interactions into account we extend Hhc by including the Fock states 
with ni — 2, 

' H^ V 
U f Hhc 



#(1+2) = 



(4) 



In the Hamiltonian (@J the coupling matrix V is proportional to J and H^\n) = W\n). Using 
ordinary perturbation theory in the parameter J/W <C 1 and projecting the eigenstates of 
back to the Hilbert space Hhc, we end up with the following effective Hamiltonian, 



H sc = — J £ ( a z+i a < 




h.c.) . (5) 



For finite W the 'soft-core bosons' (SC-bosons) Hamiltonian © approximates the low-energy 
spectrum of Bose atoms much better than the HC-bosons Hamiltonian J2J). As an example, 
Fig. n compares the low-energy spectrum of the system with the spectrum of the SC- 
bosons Hamiltonian for J/W = 0.1 where, to facilitate the comparison, we parameterised the 
creation and annihilation operators by the phase 9: a\ — * a\ exp(—iW), cii — > a; exp(i^). (The 
reason for this parameterisation will become clear in a moment.) Note that some of the level 
crossings in Fig. are actually avoided crossings. We found, however, that in the relevant 
region J/W < 0.1 the gaps of the avoided crossings can be safely neglected. This result also 
means that the single-particle quasimomentum numbers ki, i = 1, N, are still good quantum 
numbers to label the eigenstates of the Hamiltonian 10. 

3. We proceed with the Bloch dynamics. Using a gauge transformation the static forcing 
can be presented as a periodic driving of the system with the Bloch frequency wb, 

H SC {t) = H HC (t) + H mt (t) , (6) 



( 1 )The derivation of the effective Hamiltonian presented in this work is similar to that for the effective 
Hamiltonian of composed bosons in the frame of the negative-u Fermi-Hubbard model (see, Ref. [13], for 
example). 
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Fig. 1 - Spectrum of the SC-bosons Hamiltonian JSJ parameterised by the phase 8 (solid lines) as 
compared to the low-energy spectrum of the Hamiltonian @ (dashed lines). The parameters are 
L = 6, N = 3, and W = 10 J. 



H HC (t) 

Hint (t) 



in 



e ^ Bt aj +1 ai + h.c. 



W 



J 2 ( 

~ ^7 H [e- l2wBt a| +1 n/a ; _i + h.c. 



Then the dynamics of the system is conveniently described by the evolution operator over one 
Bloch period Usc{Tb), where 



Usc(t) = exp 



ft Jo 



H SC {t)dt 



(7) 



and the hat over the exponent sign denotes the time odering. For W — oo (the HC-bosons 
limit) the explicite form of the evolution operator is given by 



u HC (t) = f t B(t)f , 



(8) 



where the unitary operator T represents the transformation from the Wannier basis |n) = 
\m, . . . ,ni) to the Bloch basis |k) = \kx, . . ■ , ki) ( 2 ) and the matrix of the operator D(t) is 
diagonal in the Bloch basis, 



(k|D(i)|k) = exp 



J 

'ft Jo 



t N 



27rfcj 
L 



uiBt I dt 



(9) 



From (18191) it follows immediately that Uhc{Tb) = I and, hence, BO of HC-bosons are 
perfectly periodic. In contrast, for SC-bosons the Bloch evolution operator Usc{Tb) differs 
from the identity operator and, thus, BO are generally not periodic. 

( 2 )The elements of T are given by the determinant of the N X N matrix A with Aij = L -1 / 2 exp(ikilj /L), 
where ki and lj are the occupied orbitals. 
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Fig. 2 - Normalised mean momentum (p(t) — > p{t)/po, po = NJMd/h) of strongly interacting 
Bose atoms for weak (Fd/J — 0.1, upper panel), moderate (Fd/J = 1 middle panel), and strong 
(Fd/J = 4, lower panel) forcing. The other parameters are L = 8, N = 3, and W = 10 J. 



Addressing the dynamics of interacting Bose atoms one has to distinguish between strong, 
Fd ^> J, and weak, Fd -C J, forces. We discuss a weak force first. In this case the phase 
9 = ujBt in Eq. JfjJ) can be considered as a slowly varying parameter and, hence, the system 
adiabatically follows the instantaneous levels shown in Fig.^ Then the matrix of the evolution 
operator Usc(Tb) is diagonal in the eigen-basis of the Hamiltonian Hsc- For J/W <C 1 this 
basis practically coincides with the eigen-basis of Hhc, i-e., with the quasimomentum Fock 
states |k). Thus 

<k| U SC (T B ) | k) = expHe(k)] , (10) 

where 0(k) = j- J^ B E\^{t)dt are the phases acquired by the states during the time evolution. 
Since the energy differences between the levels of the HC-bosons and SC-bosons Hamiltonians 
are essentially given by the diagonal elements of Hi nt (t) in the Bloch basis, we finally have 

9(k) = ""lS (k| ^^ +1 " ;|k> - 

Assuming the experimentally relevant situation where the system is initially in its ground 
state, the result <|10lll[l means that BO of SC-bosons are periodic in the limit of small F. 
This regime of the system's dynamics is illustrated in the upper panel of Fig. [3 showing the 
mean atomic momentum of N = 3 atoms in L = 8 wells (periodic boundary conditions) for 
J/W = 0.1 and Fd/J = 0.1. 

When the static force is increased, the adiabatic approximation breaks down and the 
Bloch dynamics of the atoms becomes rather complicated. An example is given in the middle 
panel in Fig. [3 where the static force Fd/J =1. It is seen that the mean momentum 
p(t) undergoes a complex quasiperiodic process, involving many different frequencies. These 
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Fig. 3 - Quasienergy spectrum of the system for different magnitudes of the static force. The system 
parameters are L = 8, N = 3, and W = 10J. 



frequencies are obviously given by the quasienergies E, which are defined by the eigen-phases 
G of the evolution operator, 



U sc (T B )m = exp(-*e)|$) 



(12) 



through the relation 8 = ETs/h — 2irE /Twjb- For the considered system of N = 3 atoms 
in L = 8 wells the quasienergies E are plotted in Fig. [3] as the functions of F, the strength 
of the static force. The transition from the adiabatic to the non-adiabatic regime, reflected 
in the removed degeneracy between the quasienergy levels, is clearly seen( 3 ). It is also an 
appropriate place here to mention the particle-hole symmetry of the system. In terms of the 
evolution operator it means that the quasienergy spectum of N = 3 atoms in L = 8 wells 
coincides (up to an irrelevant global shift) with the spectrum of TV = 5 atoms in L = 8 wells. 

The further increase of F simplifies the Bloch dynamics again, as can already be concluded 
from the structure of the quasienergy spectrum for Fd ^> J (see Fig- EJ - Let us show that in 
this region BO of SC-bosons are periodically modulated with the period 



Tw 



2irhW 



W Fd 

1~T 



(13) 



(see lower panel in Fig. EJ - Indeed, for strong forcing the Bloch period Tb is the shortest time 
scale in the system and, hence, one can separate the slow (modulation) and fast (BO) dy- 
namics. Introducing the 'slow' wave-function |W(£)) , = UHc(t)\^(t))> the Schrodinger 
equation reads 



ih- 



1 

Tb 



dt 



{t)H int (t)U HC (t)dt , 



(14) 
(15) 



( 3 )Still, since the global quasimomentum is conserved, each of the quasienergy levels is L-fold degenerate. 
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Furthermore, because the quantities (n| J^, hi+ifii\n) are integers, the relation (|16|l ensures 
perfectly periodic modulations of the BO, 



The constants A and B in Eq. I|18|l depend on the filling factor n = N/L and one gets the 
most pronounced modulations for n ~ 1/2 ( 4 ). It is also worth of mentioning that the above 
analysis implicitely assumes Fd < W. Indeed, for Fd w W the static forcing resonantly 
couples the Fock states belonging to Hhc to the states with m =2 and our approach (based 
on ordinary perturbation theory) is not valid. For a particular case of the Mott-insulator 
initial state (n = 1) the dynamical responce of the system to the resonant static force was 
analysed in Refs. [15-17]. 

4- In conclusion, we have studied the effect of residual atom-atom interactions on the Bloch 
dynamics of strongly interacting Bose atoms. The characteristic energy of these interactions 
is given by the ratio of the squared hopping matrix element J and the microscopic interaction 
constant W , and is usually neglected when one discusses the equilibrium properties of the 
system (i.e., one uses the hard-core bosons model [18,19]). The effect of residual interactions, 
however, accumulates in time when it concerns the dynamics. In particular, it has been shown 
that they may modulate BO of Bose atoms. We have provided a complete description of the 
dynamical regimes of the system in terms of the Floquet-Bloch evolution operator (evolution 
operator over one Bloch period Tb = 2irh/Fd). The matrix of this operator is found to be 
diagonal in the Bloch basis in the case of a weak forcing, Fd <§; J, diagonal in the Wannier 
basis in the opposite case Fd 3> J, and non-diagonal in either of basis for a moderate forcing, 
Fd ~ J. As the result, BO of the strongly interacting Bose atoms are periodic for a weak 
static force, periodically modulated with the period TV = 2irKW/ J 2 for a strong static force, 
and involve multiple time scales for a moderate static force. 

It is interesting to estimate the modulation period TV in the typical laboratory experiment 
on BO of Bose atoms in the gravitational field. Taking, as an example, rubidium atoms in a 
3D optical lattice with d = 0.405^m {Er/2-kTi = 3558Hz), V z = 6E R , and Vj_ = 30E R , we 
have J 2 /W = 0.1 2 /T3 « 0.008£ fl , and the Stark energy Fd = 0.2AE R . Thus the modulation 
period TV = 30Tb = 35ms. This time is easily accessible in the present day experiments. It 
should be noted, however, that the results of the present work refer to a finite lattice with 
periodic boundary conditions, while in practice one typically uses a harmonic confinement of 
the system. The analysis of the Bloch dynamics of Bose atoms for these different boundary 
conditions will be the subject of a separate paper. 
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